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Abstract 

The normalization condition, average values and reduced distribution functions can 
be generalized by fractional integrals. The interpretation of the fractional analog of 
phase space as a space with noninteger dimension is discussed. A fractional (power) 
system is described by the fractional powers of coordinates and momenta. These 
systems can be considered as non-Hamiltonian systems in the usual phase space. The 
generalizations of the Bogoliubov equations are derived from the Liouville equation for 
fractional (power) systems. Using these equations, the corresponding Fokker-Planck 
equation is obtained. 
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1 Introduction 

Fractional integrals and derivatives [1] have many applications in statistical mechanics and 
kinetics [21 [3]. The generalization of the Fokker-Planck equation can be used to describe 
kinetics in the fractal media. It is known that the Fokker-Planck equation can be derived from 
the Liouville and Bogoliubov equations [51E1[7|. The Liouville equation is obtained from the 
normalization condition and from the Hamilton equations. The Bogoliubov equations can 
be derived from the Liouville equation and from the definition of the average value. In this 
paper, the generalized Fokker-Planck equation is obtained from the Liouville and Bogoliubov 
equations for fractional (power) systems. For this aim, we use fractional generalizations of 
thenormaUzation condition and the average values [HI El EDI flT] . 

In the paper, we suggest the physical interpretation of integrals of noninteger order. 
The fractional integral is considered as an integral on the fractal or noninteger-dimensional 
space. This interpretation is connected with the definition of noninteger dimension. We 
prove that fractional integration can be used to describe processes and systems on fractal. 
The physical values on fractals can be "averaged" and the distribution of the values on 
fractal can be replaced by some continuous distribution. To describe the distribution on the 
set with noninteger dimension, we use the fractional integrals. The order of the integral is 
equal to the fractal Hausdorff dimension of the set. The consistent approach to describe 
the distribution on fractal is connected with the mathematical definition of the integrals 



on fractals [THl [IH [15]. It was proved [I2] that integrals on net of fractals can be 
approximated by fractional integrals. In Ref. [El [9], [10, [H] , we proved that fractional integrals 
can be considered as integrals over the space with noninteger dimension up to a numerical 
factor. We use the well-known formulas of dimensional regularizations [16]. There is an 
interpretation that follows from the fractional measure of phase space [H |9l jTOl |TT] , which is 
used in the fractional integrals. The fractional phase space can be considered as a space that 
is described by the fractional powers of coordinates and momenta. Using this phase space, 
we can consider some of the non-Hamiltonian systems as generalized Hamiltonian systems 
[H [9l [TOl [11]. The fractional systems can be described as exitations of the fractal medium 

In Sec. 2, we consider the Hausdorff measure, the Hausdorff dimension and the integration 
on fractals to fix notation, and provide a convenient reference. The connections of the 
integration on fractals and the fractional integrals are discussed. The fractional average 
values and reduce distribution functions are defined. In Sec. 3, we derive Fokker-Planck 
equations from the Liouville equation for fractional (power) systems. A short conclusion is 
given in Sec. 4. 



2 Integration on fractal and fractional integration 
2.1 Hausdorff measure and Hausdorff dimension 

Fractals are measurable metric sets with a noninteger Hausdorff dimension. To define the 
Hausdorff measure and the Hausdorff dimension, we consider a measurable metric set {W, ^h) 
with W C M". The elements of W are denoted by x,y,z,..., and represented by n-tuples 
of real numbers x = (xi, X2, a;„) such that W is embedded in M". The set W is re- 
stricted by the following conditions: (1) W is closed; (2) W is unbounded; (3) W is regular 
(homogeneous, uniform) with its points randomly distributed. 
The diameter of a subset E G W gMJ^ is 

diam{E) = sup{G{x,y) : x,y E E}, 

where G{x^ y) is a metric function of two points x and y G W . 

Let us consider a set {Ei} of subsets Ei such that diam{Ei) < e Wi, and W C IJi^i ^i- 
Then, we define 

aEi,D) = uiD)[dtam{Ei)f. (1) 

The factor uj{D) depends on the geometry of E^. If {Ei} is the set of all (closed or open) 
balls in W, then 

The Hausdorff dimension D of a subset E G W is defined [17] by 

D = dimuiE) = sup{d e R : finiE, d) = oo} = in{{d e M : UniE, d) = 0}. (3) 
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From the definition (I3l), we obtain 

(1) ^lH{,E, d) = Ohr d> D = dimniE); 

(2) fjiniE, d) = oo ior d < D = dimu^E). 

The Hausdorff measure of a subset E dW is defined [T71 [18] by 



^h{E,D) = uj{D) lim miy^[diam{Ei)f. (4) 



diam(Ei)^0 {EA 

i=l 

Note that j2H{XE,D) = X^j^niE, D), where A > 0, and XE = {Ax, x G 



2.2 Function and integrals on fractal 

Let us consider the functions on W: 



f{x) = J2PiXEM, (5) 



i=l 



where xe is the characteristic function of E: Xe{x) = 1 if x G -E, and Xe{x) = if x ^ E. 
The Lebesgue-Stieltjes integral for ([5]) is defined by 



„ oo 

/ fd^ = y2(^if^H{Ei 



(6) 



Therefore 



/ f{x)dfiH{x)= hm y^f{xi)i{E,,D)=uj{D) lim V /(x,)[rfmm(E,)]^. (7) 

diam{Ei)—fO diam{Ei)—*Q ' ' 

Ei Ei 

It is possible to divide M" into parallelepipeds 

Ei^...i„ = {{xi, ...,Xn) eW : Xj = {ij - l)Axj + aj, < aj < Axj, j = 1, ...,n}. (8) 
Then 

n n 

di,H{x)= hm aE^.,„,^.D)= hm "Q^^^^)''^" = H 

diam(_Bi i„)-^0 diam{Ei 

The range of integration W can be parametrized by polar coordinates with r = G(x, 0) and 
angle VL. Then can be thought of as a spherically symmetric covering around a center 
at the origin. In the limit, function C,{Er^n, D) gives 

d^H{r,n)= lim ^{Er,n,D) = dn^-^r^-^dr. (10) 
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Let us consider f{x) that is symmetric with respect to some point xq G W, i.e. f{x) = 
const for all x such that G{x, Xq) = r for arbitrary values of r. Then the transformation 



W : x^x' = x-Xo (11) 



can be performed to shift the center of symmetry. Since W is not a linear space, ffTTj) need 
not be a map of W onto itself. Map ffTTl) is measure preserving. Using ffTOl) . the integral 
over a D-dimensional metric space is defined by 

2 D/2 POO 

fdf^H = :^,rj^, / f{r)r^''dr. (12) 

This integral is known in the theory of the fractional calculus [T]. The right Riemann- 
Liouville fractional integral is 



^^f^') = Y{D)j i^-^)''''fi^)d^- (13) 
Equation (fT2|) is reproduced by 

Relation f|T^ connects the integral on fractal with the integral of fractional order. This result 
permits to apply different tools of the fractional calculus [1] for the fractal medium. As a 
result, the fractional integral can be considered as an integral on fractal up to the numerical 
factor r(D/2)/[27r^/2r(D)]. 

Note that the interpretation of fractional integral is connected with the fractional dimen- 
sion [H [9l [11]. This interpretation follows from the well-known formulas for dimensional 
regularization [16] . The fractional integral can be considered as an integral in the noninteger- 
dimensional space up to the numerical factor T{D/2)/[2'k^/'^T{D)]. It was proved [H] that 
the fractal space-time approach is technically identical with the dimensional regularization. 

The integral defined in ([7j) satisfies the following properties: 

(1) Linearity: 



afi + hf2)diiH = a fidfin + b f2dfiH, (15) 
X Jx Jx 

where /i and /2 are arbitrary functions; a and b are arbitrary constants. 
(2) Translational invariance: 

f{x + Xo)dfiH{x) = / f{x)d^H{x) (16) 
x Jx 

since dfini^ — Xq) = dfini^) as a consequence of homogeneity (uniformity). 
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(3) Scaling property: 



f{ax)dfiH{x) = a ^ f{x)d^H{x) (17) 

X J X 



since d^H{x/a) = a ^d^nix). 

It has been shown [I6j that conditions (fT5l) - (fT7j) define the integral up to normalization [16]. 



2.3 Multi-variable integration on fractal 

The integral f[T^ is defined for a single variable, and not multiple variables. It is useful for 
integrating spherically symmetric functions. This integral can be generalized for the multiple 
variables by using the product spaces and product measures. 

Let us consider the measure spaces (Wk, fJ.k, D) with k = 1,2,3, and form a Cartesian 
product of the sets Wk producing the space W = Wi x W2 x W^. The definition of product 
measures and the application of the Fubini's theorem provides a measure for W as 

(/ii X /i2 X /i3)(l^) = ^^l{W^)|l2{W2)^l3{W3). (18) 

The integration over a function / on the product space is 

/(r) rf/ii X /i2 X /i3 = / / f{Xi,X2,X3) dl2i{xi)dl22{x2)dlJ,3{x3). (19) 



In this form, the single- variable measure from f|T2l) may be used for each coordinate Xk, which 
has an associated dimension ak- 

dHk{xk) = j7z:r\xk\'^''~^dxk, k = 1,2,3. (20) 

1 [ak/2) 

The total dimension of W = Wi x W2 x W3 is D = ai + a2 + 03. 

Let us reproduce the result f|T2l) from (fT9|) . We take a spherically symmetric function 
/(r) = /(xi, X2, Xs) = /(r), where = {xiY + (0:2)^ + (xa)^. Equation (fT9|) becomes 



c//ii(xi)c//i2(x2)c?/i3(x3)/(Xi, X2, X3) 

27r"i/2 27r"2/2 2Tr"^/'^ 



dr d(h rf^ J3(r,</)) r"i+"2+"3-3 



X 



r(ai/2)r(a2/2)r(«3/2) 

X (cos0)"^"^(sin0)"2+°3-2(sin0)"3-V(r), (21) 

where J3(r, 0) = sin0 is the Jacobian of the coordinate change. 
To perform the integration in spherical coordinates {r,(j),6), we use 

^^m-^xcos-.c^.^I^/ffi^, (22) 

2r([^ + z/]/2)' ^ ^ 



5 



where > 0, u > 0. Then Eq. (I2T1) becomes 

27r^/2 f 

diii{xi)dii2{.X2)dii:i{x^)f{r) = J f{r)r ^dr. (23) 

This equation describes integration over a spherically symmetric function in D-dimensional 
space and reproduces result f|T2l) . 

2.4 Probability distribution on fractal 

The probability that is distributed in the three-dimensional Euclidean space is defined by 



PsiW) = / pir)dV,, (24) 
Jw 

where p(r) is the density of probability distribution, and dV^ = dxdydz for the Cartesian 
coordinates. 

If we consider the probability that is distributed on the measurable metric set W with 
the fractional Hausdorff dimension then the probability is defined by the integral 

Pd{W) = [ p{r)dVD, (25) 
Jw 

where D = diniHiW) = ai + 02 + «3, and 

dVo = dfii{xi)dn2{x2)dfi3{x3) = €3(0, r)dV3, , (26) 

There are many different definitions of fractional integrals pj. For the Riemann-Liouville 
fractional integral, the function 03(0, r) is 

where x, y, z are the Cartesian coordinates, and D = ai + 02 + 0:3, < D < 3. As the 
result, we obtain the Riemann-Liouville fractional integral jlj in Eq. fl25l) up to numerical 
factor 87r^/2_ Therefore, Eq. can be considered fractional generalization of Eq. 

[])• 

For p(r) = p(|r|), we can use the Riesz definition of the fractional integrals [1]. Then 

C3(Ar) = ^5^^^§^kr-, (29) 
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Note that 



Therefore, we suggest to use 



Jmi_C3(£>,r) = (4i»/^)-'. (30) 



C3(Ar)^H9|f^|r|- (31) 

Definition (13T]) allows us to derive the usual integral in the limit D ^ (3 — 0). 

For D = 2, Eq. (!25|) gives the fractal probability distribution in the volume. In general, 
it is not equivalent to the distribution on the two-dimensional surface. Equation (128|) is 
equal (up to numerical factor Svr^/^) to the integral on the measurable metric set W with 
Hausdorff dimension dimH{W) = D. To have the usual dimensions of the physical values, 
we can use vector r, and coordinates x, j/, z as dimensionless values. 



2.5 Fractional average values 

To derive the fractional analog of the average value, we consider the fractional integral 
for function /(x). If function /(x) is equal to the distribution function p(x), then we can 
derive the normalization condition. If the function f{x) is equal to the multiplication of 
distribution function p(x) and classical observable A{x), then we have the definition of the 
fractional average value. 

The fractional generalization of the average value [SI [HI [IHl E] can be present by 

< A >,= {IlAp){y) + {I^Ap){y), (32) 

where 

+^ T{a)J_^{y-xy~-' T{a)J^ {x - yY^- ^ > 

For a = 1, Eq. (!32|) gives the usual average value. 

The fractional average value (l32l) can be written [H [9l [TOl [11] as 

1 

<A>„=-/ pp)(y-x) + (Ap)(y + x)]rf/i„(x), (34) 

J -oo 

where 

apa[x) — -p/ N — -p/ \- I'J'-'J 

Here, we use 

x" = S(?n(x)|x|", (36) 

where the function sgn{x) is equal to +1 for a; > 0, and —1 for a; < 0. 

Let us introduce notations to consider the fractional average value for phase space. 
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(1) The operator T^j, is defined by 

T^J{...,Xk,...) = ^ (^/(..., 4 -Xfc, ...) + /(..., x'fc + Xfc,...)). (37) 

For /c-particle, which is described by coordinates qks and momenta pks, (s = l,...,m), the 
operator T[k] is 

m=T,,,T,^,...T,,^T,^^. (38) 

For the n-particle system, we define the operator T[l, ...,n] = T[l]...T[n]. 

(2) The operator I"^ is defined by 

/ + 00 
f{xk)d^^c.{xk). (39) 
-oo 

Then the fractional integral (13^ can be rewritten in the form 

<A>^=i:T,A{x)pix). 
The integral operator /"[A;] = IgJp,,-Ig,Jp,^ is 

I'^[k]f{(ik,Pk) = j f{cik,Pk)dfXa{qk,Pk), (40) 



where 



c?/U„(qfc, pfc) = {aT{a)) '"^dq'^^ A c/p^^ A ... A dq"^^ A dp; 



km' 



For the n-particle system, we use J"[l, n] = /"[1]...J"[?t,]. 

The fractional average values for the n-particle system is defined [H [9l [TOl |1T] by 

< A >,= r [1, n]T[l, n]Ap„. (41) 

In the simple case (n = m, = 1), the fractional average value is 

/ + 00 P + OO 
/ c//i«(g,p) TgTpA{q,p)p{q,p). (42) 
-oo ^ — oo 

Note that the fractional normalization condition [8],[9],[T0l[lT] is a special case of this definition 
of average values < 1 >«= 1. 
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3 Fokker-Planck equation from Liouville equation 

Let us consider a system with n identical particles and the Brownian particle. The distribu- 
tion function of this system is denoted by /0„+i(q, p, Q, -P, t), where 

q = (qi7 •••) qn), Q/fc = (O'fcl, '?fcm), 
p = (Pl, ...,p„), Pfc = (Pfcl, •••,Pfcm) 

are coordinates and momenta of the particles; Q = {Qs) and P = (Pg) {s = l,..,m) are 
coordinates and momenta of Brownian particles. The fractional normalization condition 
P El [ini [n] has the form 

r[l,...,n,n + l]p„+i = 1, (43) 

where 

pn+i = T[l, ...,71,71 + l]p„+i(q, p, Q, P, t). (44) 
The reduced distribution function for the Brownian particle is 

ps(Q, P, t) = r [1, n]p„+i(q, p, Q, P, t). (45) 

The distribution p„+i satisfies the Liouville equation [81 [Q], \T0\ [TT]: 

i{Ln + LB)Pn+l = 0, (46) 



dpn+1 



dt 

where L„ and Lb are Liouville operators such that 



(47) 



The forces and /s, and the velocities and are defined by the Hamilton equations 
of motion. For kth particle, 

^ = G^(q,p), ^ = P,^(q,p,Q,P), (49) 
where k = 1, ...,n. The Hamilton equations for the Brownian particle are 

-|^=^?,(g,P), ^ = /,(q,p,Q,P). (50) 
For simplification, we suppose 

G]=plJm, gs = P:/M, (51) 
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where M ^ m. 

Let us use the boundary condition in the form 



hm p„+i(q,p,Q,P,t) = p„(q,p,Q,T)pB(Q,P,t), 

t— +— oo 



where pn is a canonical Gibbs distribution function 

Pn(q, P, Q, T) = exp p J^- H{<i, p, Q) 
Here, i?(q, p, Q) is a Hamilton function 

H{q, p, g) = i7„(q, p) + ^ Ueicik, Q) 

k=l 



(52) 



(53) 



(54) 



where Hn is a Hamiltonian of the n-particle system, and f/^ is an energy of interaction 
between particles and the Brownian particle. If we use Eqs. fH9|) and flSTj) . then 



^n(q,p) = + 

fc,s k<l 



(55) 



The boundary condition (1521) can be realized [12] by the infinitesimal source term in the 
Liouville equation 

^^ti _ i(^L^ + LB)Pn+l = -e{pn+l - PnpB)- (56) 

Integrating ( l56l) by /"[I, J^], we obtain the equation for the Brownian particle distribution 



dt ^ dQ 

s=l 

The formal solution [TSl has the form 



dpB , ^ djgspB) , far, ^ ^(/.P. 



'n+lj 



9Pf 



0. 



(57) 



Pn+l 



^0 

{t)=e dT e""e-'"(^"+^^)pB(t + r)p, 

^ — oo 



or 



P„+i(t) = PB{t)pn - j e^^e-^^(^"+^^) - «(L„ + L^)^ PB(t + r)p„. 
Substituting ( l59l) into ( 1571) . we get 
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(58) 



(59) 



5PB , d{gspB) , OpB far. 1/ / ~ ^ 

:^ [1, (jsPn)- 



dt 



s=l 



s=l * 



Q /-O /a \ 



Note that /"[I, ■■,n\fspn can be considered as an average value of the force fg- This average 
value for the canonical Gibbs distribution (l53l) is equal to zero. The last term can be 
simplified. Using 

'./i^Vn, (61) 



ag^ kT' 



where f'f^ is a potential force 



/.'" - -ig. (62) 



we get 



It can be proved by interactions that the term 



dt dQ 

in the integral of ( !60|) does not contribute. Then 



dpB ^ djgspB) ^gg^ 



dps 
dt 



This equation is a closed integro-differential equation for the reduced distribution function 
Pb- Note that force fs can be presented in the form 

/. = /i^^ + /i"\ 

where /i^^ is a potential force (1^2]) . and /i"^ is a non-potential force that acts on the Brownian 
particle. For the equilibrium approximation P ~ (MkTY^'^, iL^ ~ M~^/^, iL„ ~ m~^/^, 
and M ^ m, we can use perturbation theory. 

Using the approximation pB{t + r) = pB{t) for Eq. fl64|) . we obtain the Fokker-Planck 
equation for fractional power systems 

-df + l^ ^Qf- + 1^0^ [j + Iss'Ps'PBit) J = 0, (65) 

S = l " S = l i \ ' A / 
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where 

-0 



7,V = /3Mr[l,...,n] / dr e^Vse-^^^"/.'Pn. (66) 

J —oo 

7^, = /3MJ-[1, n] f dr e^/.e-^^^^/if^p^. (67) 

J — oo 

If/, = /i^\then 7L'=7L'- 

Let us consider the one-dimensional stationary Fokker-Planck equation fl65p with 

d{9spB)/dQ: = 0. 

Then 



Obviously, we get the relation 

Md{^'{P)pBit)) 



13 

Assuming that the constant is equal to zero, we get 



+ -f^{P)PpB{t) = const. (69) 



9P° M 
or, in an equivalent form 

(9P° 7i (P)M' 

The solution is 

„2 



(71) 



ln[^\P)PB{t)] = I ^^^^dP- + const. (72) 



As the result, we obtain 



where the coefficient is defined by the normalization condition. Equation (I73|l describes 
the solution of the stationary Fokker-Planck equation for the fractional (power) system. The 
special cases of f l73|) can be derived as done in Ref. 
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4 Conclusion 



In this paper, the fractional generahzations of the average value and the reduced distribution 
functions are used. The generalization of the Liouville and Bogoliubov equations are derived 
P, El Uni [H] from the fractional normalization condition. Using these equations, we obtain 
the Fokker-Planck equation for fractional (power) systems. 

The Liouville, Bogoliubov and Vlasov equations for fractional systems [H [9l [TOl [11] can 
be considered as equations in the noninteger-dimensional phase space. For example, the 
systems that live on some fractals can be described by these equations. Note that the frac- 
tional systems can be presented as special non-Hamiltonian systems [SI El HDl E] . Fractional 
oscillators can be interpreted as elementary excitations of some fractal medium with non- 
integer mass dimension [SI El [lOl E]- The fractional (power) systems are connected with 
the non-Gaussian statistics. Classical dissipative and non-Hamiltonian systems can have 
the canonical Gibbs distribution as a solution of the stationary equations [20l [21]. Using 
the methods [201 [21] , it is easy to prove that some fractional dissipative systems can have 
the fractional analog of the Gibbs distribution (non-Gaussian statistic) as a solution of the 
stationary equations for fractional systems. 
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